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This  report  is  carposed  of  two  parts  covering  the  response  of  spinal 
ligaments  to  cyclic  creep  and  a  mathematical  model  of  those  ligaments  to  a 
specified  strain  history.  This  concludes  a  three-year  endeavor  to  deter¬ 
mine  the  properties  of  spinal  ligaments  frcm  rhesus  monkey,  baboon,  chimpan¬ 
zee  and  human  specimens. 

>  Four  ligaments  of  the  spine  of  rhesus  monkey  and  baboon  were  selected 
for  testing  in  Part  1'  with  test  sanples  taken  at  different  vertebral 
levels.  The  cyclic  creep  tests  on  ligaments  from  three  rhesus  monkeys  and 
two  baboons  were  preceded  by  single  extension  failure  tests  on  ligaments 
from  two  rhesus  monkeys  to  establish  a  range  of  static  failure. 

As  described  -in  Part  2,  a  mathematical  constitutive  equation  was 
developed  using  a  linear  viscoelastic  model  for  individual  fibers  and  a  dis¬ 
tribution  function  to  include  the  effects  of  fiber  distribution,  orienta¬ 
tion  and  initial  length.  Material  constants  and  the  distribution  function 
were  determined  by  numerical  solution  of  the  constitutive  integral  equation 
using  the  fastest  strain  rate  data.  Analytical  and  experimental  results 
were  compared  for  selected  human  spinal  ligament  sanples. 

This  research  completes  an  initial  investigation  of  the  mechanical 
properties  of  spinal  ligaments  and  ligaments  from  the  lower  extremities. 
Research  of  this  nature  is  necessary  to  understand  the  mechanism  of  injur¬ 
ies  in  airnrewmen  that  result  from  exposure  to  vibratory  loading  environ¬ 
ments  and  high  stress  or  strain  rates  during  escape  and  crash  episodes. r 
Determining  the  tissue  response  differences  between  different  primatel 
species  will  aid  in  the  selection  of  animal  models  and  in  development  of 
the  interspecies  scaling  techniques  to  establish  criteria  and  of  preventive 
measures  for  the  safety  of  aircrews. 
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Biological  tissues  are  viscoelastic  so  that  their  responses  to  a  cur¬ 
rent  loading  or  deformation  are  affected  by  their  mechanical  environment 
before  the  current  time.  Loading  and  deformation  of  tissues  occur  con¬ 
tinuously  by  exercise  or  inputs  fran  outside  the  body.  The  continuous 
nature  of  demands  on  tissues  and  the  dependence  of  tissue  response  on  pre¬ 
vious  demands  lead  to  significant  questions  concerning  the  effects  of 
repeated  loading  on  tissue  deformation  and  failure.  For  exaxtple,  what  are 
the  effects  of  vigorous  exertion  or  severe  mechanical  inputs  on  the  poten¬ 
tial  for  injury  of  a  pilot  in  an  emergency  situation?  In  studying  and 
modelling  the  biodynamic  response  of  the  pilot,  hew  are  the  responses  of 
tissues,  and  the  modelling  elements  which  represent  them,  affected  by 
repeated  loading? 


The  load  relaxation  response  to  a  few  extension  cycles  has  been  used 
to  evaluate  the  validity  of  constitutive  models  for  soft  connective,  tis¬ 
sues;  typical  among  such  studies  are  the  works  of  Haut  and  Little1  and 
Woo,  .Gamez,  and  Akeson  which  are  based  on  the  viscoelastic  theories  of 
Fung.  Lanir  has  proposed  a  model  for  soft  connective  tissues  that 
incorporates  the  approach  of  Fung1  with  a  distribution  of  tissue  fiber 
lengths.  A  similar  model  is  presented  and  discussed  in  Part  2  of  this 
report.  Lanir4  predicts  load  relaxation  responses  to  constant  and  cyclic 
extensions  and  the  creep  extension  responses  to  constant  and  cyclic  load¬ 
ing;  however,  no  experimental  evaluations  of  tissue  behaviors  have  been 
made  of  Lanir' s  model,  nor  has  the  model  presented  in  this  report  been  com¬ 
pared  with  experimental  data  for  cyclic  creep. 


-The  creep  extension  of  human  digital  tendon  was  studied  by  Cohen  et 
al. ,  with  particular  attention  to  the  activation  energy  for  the  visco¬ 
elastic  mechanism.  They  found  creep  strain  rate  to  be  roughly  proportional 
to  stress  raised  to  a  positive  power  less  than  one.  The  creep  extension  of 
passive  cardiac  muscle  was  measured  by  Pinto  and  Fung  and  was  found  to 
be  roughly  a  linear  function  of  logarithm  of  time.  The  bending  deformation 
and  recovery  due  to  application  and  removal  of  loading  on  fresh  human  lum¬ 
bar  spines  were  recorded  by  Tvramey  and  Taylor. ' 


The  previous  research  on  creep  of  connective  tissues  was  very  limited 
in  extent  and  included  no  data  on  creep  deformation  due  to  cyclic  loading 
or  on  the  effects  of  repeated  loading  on  tissue  failure. 


The  current  study  reported  in  Part  1  was  undertaken  to  obtain  initial 
data  for  the  effects  of  repeated  loading  on  the  deformation  and  failure  of 
spinal  ligaments.  The  approach  to  this  research  was  modelled  after  the 
study  of  metal  fatigue  in  vhich  material  samples  are  failed  by  single  exten¬ 
sions  or  by  repeated  loading,  and  the  failure  loads  are  normalized  by  sam¬ 
ple  area  to  obtain  plots  of  stress  versus  number  of  cycles  to  failure. 
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The  goals  of  Part  1  of  the  research  were  to: 

a.  Study  the  effects  of  cyclic  loading  on  the  failure  of  spinal 

ligaments; 

b.  Study  the  sattple  extensions  due  to  the  cyclic  loading,  i.e., 

cyclic  creep. 

The  goal  of  Part  2  was  to  model  selected  spinal  ligaments  in  terms  of 
their  anatomical  microstructure  (including  fibral  composition  and  orien¬ 
tation)  as  well  as  their  biomechanical  response. 


PART  1 


CYCLIC  LOADING  OF  SPINAL  LIGAMENTS  FROM  RHESUS  MONKEYS  AND  BABOONS 


METHODS  AND  MATERIALS 


Spinal  ligament  samples  from  two  rhesus  monkeys  were  pulled  to  failure 
in  single  extensions,  and  spinal  ligament  samples  from  three  additional  rhe¬ 
sus  monkeys  and  two  baboons  were  subjected  to  cyclic  loading  and  failure. 
The  methods  used  in  preparing  the  spinal  ligaments  for  testing  and  in  sam¬ 
ple  gripping  were  the  same  as  in  previous  mechanical  testing  of  spinal  liga¬ 
ments.  That  is,  the  spinal  column  was  removed  from  an  eviscerated  pri¬ 
mate,  sectioned  into  vertebral  pairs  for  testing  and  histology  according  to 
Figure  1,  and  sectioned  further  into  bone- ligament-bone  sanples  containing 
either  the  anterior  longitudinal  ligament  (A.L.L. ),  posterior  longitudinal 
ligament  (P.L.L.),  ligamentum  flwum  (L.F.),  or  supraspinous  ligament 
(S.S.L. ).  As  in  previous  studies,  the  sanples  were  gripped  by  clamping 
the  vertebral  bone  segments  between  fixtures  with  waterproof  sanding  mesh 
of  80  grit;  the  ligament  cross-sectional  areas  were  measured  using  histo¬ 
logical  sections. 

The  spinal  ligament  samples  from  two  rhesus  monkeys  were  extended  to 
failure.  These  failure  loads  were  normalized  by  ligament  cross-sectional 
area  resulting  in  failure  stress  and  were  used  to  estimate  the  failure 
loads  for  the  other  ligament  samples  that  were  loaded  repeatedly. 

Ligament  lengths  ware  approximated  as  the  bone-to-bone  distances  when 
samples  were  mounted  in  the  testing  grips  and  extended  to  onset  of  loading. 
Based  on  these  estimates  of  ligament  length,  the  deflection  rates  for  the 
single-extension  failure  tests  were  selected  to  produce  1%/s  strain  rate 
and  resulted  in  an  average  stress  rate  of  2  MPa/s  based  on  failure  stress 
and  test  duration.  Since  the  samples  stiffen,  the  stress  rate  at  failure 
would  be  higher  than  this  average. 

A  testing  protocol  was  developed  to  subject  the  ligament  sanples  to 
repeated  cycles  of  load  with  sample  loading  controlled  to  vary  between  mini¬ 
mum  and  maximum  values  that  were  less  than  the  estimated  single-extension 
failure  load.  The  stress  levels  in  this  cyclic  testing  could  only  be  esti¬ 
mated;  exact  stress  values  were  determined  after  histological  studies  of 
ligament  areas  were  conducted.  The  cyclic  tests  were  run  with  the  load 
sinusoidally  ranging  from  the  maximum  load  to  an  average  of  0.65  of  that 
maximum;  that  is  an  average  of  about  20%  cyclic  variation  in  load  fran  a 
constant  mean  load.  The  average  cyclic  stress  rate  was  22  MPa/s,  which  was 
super inposed  on  the  mean  stress  level.  The  cyclic  stress  rate  was  about 
ten  times  greater  than  the  stress  rate  in  the  failure  tests  by  single  exten¬ 
sion.  The  fact  that  the  cyclic  stress  was  superimposed  on  a  constant 
stress  makes  the  cortparison  of  stress  rates  between  test  types  difficult. 
However,  increased  loading  rate  is  kncwn  to  increase  failure  loads  in  liga- 
mentSpand  increase  the  likelihood  of  bone  failure  in  bone-ligament-bone  sam¬ 
ples. 
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Samples  from  the  first  animal  tested  were  loaded  cyclically  to  failure 
or  it  was  apparent  that  they  would  not  fail  at  that  loading  level,  and  the 
test  was  stopped.  Later,  sanples  were  cycled  at  chosen  load  levels  until 
they  failed  or  for  about  4,000  cycles  after  which  they  were  extended  to 
failure. 

The  testing  was  performed  with  a  servohydraulic  materials  testing  sys¬ 
tem,  and  the  entire  cyclic  testing  protocol  was  controlled  and  recorded  by 
digital  computer.  Once  the  sanple  was  mounted  in  the  testing  grips  and 
appropriate  testing  parameters,  such  as  sample  identification  and  loading 
limits,  had  been  entered  into  the  computer ,  the  sairple  was  extended  under 
stroke  control  to  the  minimum  load  for  cyclic  testing;  and  then  the  testing 
machine  was  transferred  to  load  control  for  cyclic  loading.  All  these 
tests  were  run  at  a  frequency  of  one  cycle  per  second.  The  tests  were 
started  in  stroke  control  to  avoid  possible  servocontrol  instability  in 
load  control  with  a  saiple  of  low  initial  stiffness. 

During  the  cyclic  testing,  load  and  extension  were  recorded  with  the 
computer  at  50  data  pairs  per  second  for  about  5  seconds  in  repeated  fif¬ 
teen  second  intervals.  The  sanpled  data  pairs  were  filtered  by  a  running 
average  of  three  points  and  scanned  for  maximum  and  minimum  values  of  load 
and  extension.  These  values,  with  corresponding  times,  were  then  written 
to  the  conputer  disc  for  storage  and  subsequent  analysis.  Thus,  the  data 
stored  in  the  computer  from  each  test  were  maximum  and  minimum  values  of 
load  and  extension,  sampled  periodically  rather  than  continuous  records  of 
these  parameters,  which  would  have  been  far  too  extensive  for  reasonable 
data  storage.  In  addition  to  the  data  stored  throughout  the  test  by  the 
computer,  the  initial  30  seconds  of  cyclic  load  and  extension  were  recorded 
at  50  data  pairs  per  second,  using  a  digital  oscilloscope  with  magnetic 
storage  discs. 

The  particulars  concerning  the  animals  and  respective  tests  are 
presented  in  Table  1. 

Table  1.  Animal  Identification  Number,  Species, 

Sex,  Weight,  and  Test  Type 


Identification 

Species 

Sex 

Weiqht,  kq 

Test  'rype 

U70 

rhesus 

male 

9.0 

failure 

U88 

rhesus 

male 

13.9 

failure 

N75 

rhesus 

female 

7.0 

cyclic 

A182 

rhesus 

male 

10.0 

cyclic 

A322 

rhesus 

male 

8.5 

cyclic 

H-40 

baboon 

male 

14.1 

cyclic 

H-54 

baboon 

male 

32.7 

cyclic 

RESULTS 


The  results  frcm  failure  testing  of  spinal  ligaments  from  two  rhesus 
monkeys  are  given  in  Tables  2  and  3.  By  knowing  the  ligament  areas  from 
histological  sections  adjacent  to  the  tested  ligaments,  their  areas  were 
interpolated  and  the  failure  loads  were  divided  by  the  interpolated  cross- 
sectional  areas  to  calculate  failure  stress.  The  failure  loads  and  liga¬ 
ment  areas  for  the  larger  rhesus  (U88)  were  generally  greater  than  for  the 
smaller  animal  (U70),  so  that  the  failure  stresses  were  similar  for  same 
ligament  and  level  in  each  animal.  Large  variability  in  the  failure 
stresses  for  supraspinous  and  posterior  longitudinal  ligaments  (S.S.L.  and 
P.L.L.,  respectively)  was  due  to  the  small  size  and  indistinct  nature  of 
these  ligaments  in  corrparison  to  ligamentum  flavum  (L.F. )  and  anterior 
longitudinal  ligaments  (A.L.L. ).  The  extensions  at  failure  for  the  P.L.L. 
were  generally  lower  than  the  other  ligaments.  The  L.F.  and  A.L.L.  sup¬ 
ported  several  times  more  load  at  failure  than  P.L.L.  or  S.S.L.  The  fail¬ 
ure  stress  and  extension  data  from  single  extensions  to  failure  will  be  sum¬ 
marized  and  coupared  to  such  data  from  cyclic  loading  later  in  this  part  of 
the  report. 

After  the  failure  tests  by  single  extension,  spinal  ligament  samples 
from  three  rhesus  monkeys  were  subjected  to  cyclic  loading.  Stress  levels 
were  chosen  which  were  below  the  stresses  that  had  caused  failure  of  corres¬ 
ponding  ligaments  in  rhesus  monkeys  U70  and  U88.  For  samples  from  the 
first  such  animal  (rhesus  N75),  the  intention  was  to  continue  the  test 
until  each  sanple  failed  during  load  cycles  at  less  than  the  single-exten¬ 
sion  failure  stress.  As  can  be  seen  from  Table  4,  seme  of  these  samples 
failed  immediately  or  nearly  so;  some  failed  after  a  few  hundred  cycles; 
and  some  endured  several  thousand  cycles.  During  the  testing  of  the  sam¬ 
ples  from  the  next  animal  (rhesus  A182),  it  was  decided  to  limit  the  number 
of  cycles  to  4,000;  and  if  the  samples  had  not  failed  by  then,  they  would 
be  pulled  to  failure  by  a  single  extension.  Thus,  for  some  of  the  liga¬ 
ments  in  the  last  two  rhesus  monkeys  (A182  and  A322),  the  failure  stresses 
were  greater  than  the  maximum  cyclic  stresses  (Tables  5  and  6). 

Cyclic  loading  was  also  performed  with  spinal  ligaments  from  two 
baboons.  The  maximum  cyclic  load  levels  were  selected  based  on  failure 
loads  from  the  rhesus  monkeys  and  were  scaled  to  the  baboons  by  considering 
the  ligament  areas.  The  cyclic  tests  with  the  baboon  ligaments  were  run 
similarly  to  those  of  the  last  two  rhesus  monkeys;  that  is,  the  samples 
were  run  for  up  to  4,000  cycles,  and  if  they  had  not  failed,  they  were 
pulled  in  a  single  extension  to  failure.  The  results  of  these  baboon  tests 
are  given  in  Tables  7  and  8.  In  Tables  4  through  8,  a  value  of  failure 
stress  is  given  unless  the  sairple  did  not  fail  during  cyclic  loading  and 
was  not  subsequently  extended  to  failure. 

The  data  from  Tables  2  through  8  are  presented  in  Figures  2  through  5 
as  "stress  versus  number -of -cycles"  plots,  commonly  used  in  studies  of 
metal  fatigue.  The  data  are  grouped  by  ligament  type,  i.e.  S.L.L.,  L.F., 
P.L.L.,  and  A.L.L. ,  and  by  location,  i.e.  thoracic  (T3-T4,  T6-T7,  and 
T9-T10)  or  lumbar  (T12-L1,  L2-L3,  and  L5-L6).  The  values  of  failure 
stresses  for  single  extensions,  either  by  design  for  the  first  two  animals 
tested  or  by  chance  due  to  scatter  in  strength  of  the  ligaments  from  the 
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subsequent  tests,  are  plotted  with  maximum  and  minimum  cyclic  stresses  and 
failure  stresses  at  the  corresponding  number  of  cycles.  The  scatter  in 
strength  was  apparent  from  these  Figures.  The  ligaments  appeared  to  be  as 
strong  after  hundreds  or  thousands  of  load  cycles  as  t  av  were  for  a  single 
extension.  That  is,  there  did  not  appear  to  be  any  significant  effect  of 
cyclic  loading  on  the  strength  of  the  spinal  ligaments  tested.  However, 
the  failure  loads  in  the  cyclic  tests,  in  comparison  to  the  single  exten¬ 
sion  tests,  may  have  been  elevated  by  the  average  cyclic  stress  rate  of  22 
MPa/s  versus  2  MPa/s  for  the  single  extension  tests. 

The  failure  stresses  are  suntnarized  in  Table  9,  and  the  deflections  at 
which  the  samples  failed  are  summarized  in  Table  10,  respectively,  for  the 
single  extensions  to  failure,  cyclic  loading,  and  these  two  test  types  com¬ 
bined.  The  samples  have  been  divided  by  region  into  thoracic  and  lumbar. 
The  variability  is  either  the  standard  deviation  for  numbers  of  samples 
that  are  four  or  greater  or  range  for  three  samples  or  less. 

In  Table  9  there  were  no  significant  differences  between  failure 
stress  levels  for  the  same  ligament  type  due  either  to  location  (thoracic 
versus  lumbar)  or  best  type  (single  extensive  versus  cyclic  loading)  with 
the  following  exceptions.  Lumbar  A.L.L.  appeared  to  have  a  higher  failure 
stress  in  single  extension  and  cyclic  loading  than  thoracic  A.L.L.  There 
were  differences  in  strength  between  ligament  types  with  the  L.F. ,  which  is 
largely  elastin,  able  to  resist  about  half  the  stress  of  the  other  liga¬ 
ments,  which  are  primarily  collagen. 

The  extensions  at  failure  also  indicated  no  significant  effect  of 
cyclic  loading  compared  to  single  extension  or  effect  of  location  (thoracic 
versus  lumbar);  i.e.,  any  such  effects  were  masked  by  variability  due  to 
other  factors.  The  mean  values  of  failure  extension  for  P.L.L.  were  con¬ 
sistently  less  than  for  the  other  ligaments. 

The  cyclic  testing  of  spinal  ligaments  between  maximum  and  minimum 
load  values  resulted  in  cyclic  extensions  which  increased  from  cycle  to 
cycle.  This  extension  is  cyclic  creep.  The  cyclic  extensions  were 
recorded  in  short-term  with  the  digital  oscilloscope  for  about  the  first  30 
cycles  and  in  long-term  with  the  digital  computer  every  15  seconds  from 
test  initiation  until  test  completion.  After  trying  several  mathematical 
functions,  these  short-  and  long-term  data  were  found  to  be  best  described 
by  a  power  function  of  time: 


where:  Sn  is  the  maximum  or  minimum  extension  normalized  to 
the  respective  value  from  the  first  cycle; 
t  is  time  measured  in  seconds; 

S  is  an  exponent. 


10 


This  function  is  a  straight  line  through  the  origin  of  a  plot  of  loga¬ 
rithm  of  normalized  extension  versus  logarithm  of  time  with  the  exponent 
(8)  as  the  slope.  By  taking  the  logarithm  of  both  sides  of  equation  1, 
the  following  equation  is  obtained  which  is  used  for  linear  regression  with 
both  short-  and  long-term  data  sets,  i.e. ,  frcm  the  digital  oscilloscope 
and  frcm  the  conputer,  respectively: 

In  Sn  =  8[ln( t) ]  (2) 

The  maximum  and  minimum  extensions  were  normalized  by  their  respective 
values  in  the  first  load  cycle  and  fit  to  equation  2  with  8  equal  to  Bl 
for  the  maximum  extensions  and  8  equal  to  B2  for  the  minimum  extensions. 

The  results  of  the  exponential  fits  are  given  in  Table  11.  The  corre¬ 
lation  coefficients  for  these  regressions  were  consistently  high,  typically 
above  0.97,  so  there  is  high  confidence  that  the  cyclic  creep  data  are  well 
fit.  The  exponents  Bl  and  B2  were  generally  greater  for  the  short-term 
creep  response  than  the  long-term,  indicating  an  initial  creep  that  was 
more  rapid  than  the  responses  over  hundreds  or  thousands  of  load  cycles. 
Although  there  was  variability  in  the  exponents,  it  was  apparent  that  the 
creep  in  the  L.F.  was  less  than  the  other  ligaments.  This  is  consistent 
with  the  fact  that  the  L.F.  is  largely  elastin,  whereas  the  other  ligaments 
are  largely  collagen. 

The  exponent  for  the  extension  associated  with  minimum  cyclic  load, 
B2,  was  generally  greater  than  the  exponent  for  the  peak  extension,  Bl. 
This  indicated  that  the  samples  were  stiffening,  with  the  difference 
between  minimum  and  peak  extension  decreasing  with  repeated  cycles  of  load. 
The  ratio  B2/B1  was  proportional  to  the  rate  of  stiffening  within  the  range 
of  numbers  of  cycles  in  this  study. 


PART  2 


MATHEMATICAL  MODELING 


The  spinal  ligaments  have  been  modeled  to  reflect  both  their  anatomi¬ 
cal  microstructure  and  fibral  ccnposition.  These  ligaments  are  conposed  of 
fibers  which  vary  in  initial  length  and  in  orientation  with  the  spinal 
column  axis.  These  variations  contribute  to  the  observed  non-linear 
response  of  the  ligaments  as  the  different  fibers  straighten  and  become 
aligned  with  the  loading  axis  at  different  strains  in  the  loading  history. 
This  observed  non-linearity  was  introduced  into  the  mathematical  model  by 
use  of  a  conplex  distribution  function  similar  to  that  suggested  by 
Lanir.  This  distribution  function  includes  the  effects  of  fiber  orien¬ 
tation,  initial  waviness,  and  distribution. 

Each  fiber  when  straightened  was  considered  to  be  a  linearly  viscoelas¬ 
tic  element  satisfying  a  constitutive  equation  of  the  form: 


(T^t) 


ln(t  -  X  +  D] 


(3) 


where  ov  is  the  fiber  stress,  i  is  the  fiber  index,  E  and  p  are  fiber 
material1  constants,  e  .  is  the  fiber  strain,  t  is  a  particular  time,  and  X, 
is  the  variable  time  for  integration. 

Consider  a  single  fiber  at  an  angle  0^  to  the  loading  axis  as  shown 
in  Figure  6.  The  fiber  is  considered  to  be  wavy  with  an  initial  length  of 

*i- 


This  initial  length  may  be  written 


=  x  (1  +  £i) 
1  cos  0^ 


(4) 


where  £.  is  the  ligament  strain  at  which  this  ith  fiber  assumes  load 
corresponding  to  a  ligament  deformation  of  Ax..  The  fiber  deformation 
becomes 


A .  =  (Ax  -Ax.)  cos  9  . 


(5) 


and  the  fiber  strain  is 


£i  = 


(6  -  cos2  0^ 
(1  + 


where  £  represents  the  observed  ligament  strain. 

Substituting  equation  6  into  equation  3  yields 


<Ti  -  E  Si  El 


If  cos  V  . 

+  p  ln(t  -  t  +  1)]  H - -  dt* 

dt(l  +  £i) 


The  lower  limit  of  integration  corresponds  with  the  time  the  ith  fiber 
assumes  load. 

The  axial  force  in  the  ligament  is  equal  to  the  sum  of  the  fiber  load 
components. 


F  =  X  <T  j_  \  cos  0  i 


where  A{  is  the  area  of  the  ith  fiber. 

A  carp  lex  distribution  function  <£(£)  will  be  introduced  to  include  the 
effects  of  angle  and  initial  length.  The  limits  of  integration  for  each 
fiber  may  be  ccnpensated  for  by  use  of  the  Heavyside  step  function 
H(  £  ) .  Substituting  equation  7  into  equation  8  and  integrating  over  the 

area  yields 

F  =  EA  f  <*(*))  S  U  +  H  ln(t  -  *£+  D]  H<£ -Y»  )  dtdv) 

0  1  0  dL  L  (9) 

where  A  is  the  total  area  of  the  ligament,  and  £  is  the  observed  ligament 
strain.  The  distribution  function  «C(£)  includes  the  orientation  angle 
6^,  preload  fiber  strain  £.  and  fiber  area  or  number  A..  The 
contribution  of  each  of  these  ractors  to  the  constitutive  re'ation  cannot 
be  independently  determined.  The  distribution  function  is  obtained  for 
each  ligament  by  solving  the  follcwing  constitutive  equation  using  the 
highest  constant  strain  rate  data. 


r£(t)  rt 

<T=E  J  <*.(«)  J  [1  +  m  ln(t  -  7:+  1)]  H(£  -VI  ) 
0  0 


at 


dX  dT^ 


(10) 
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The  distribution  function  can  be  expanded  in  a  power  series  in  the 

form 


*“l>  -  Z  vf 

m  v 


(id 


Substituting  this  expansion  into  equation  10  yields 
_  r£(t)  rt 

(T=  E  Y.  cm  J  J  [1  +  M  ln(t -'t  +  1)]  H(£ -v^)  l^dtdv^  (12) 

This  form  of  the  constitutive  equation  may  be  integrated  for  simple  strain 
histories. 

Constant  strain  rate  loading: 

A  constant  strain  rate  loading  has  a  strain  history  of  the  form 

£(t)  =  at  0  <  t  <  t,  (13) 

Substituting  into  equation  12  yields 

(T(t)  =  E  ]T  C.  f  1f)m  /  [1  +  \x  ln(t  -TL+  1)]  H(t,-^)  a  dtdr^  (14) 

m  0  ^  0 


Integrating  gives 


QT(t)  =  E  X 

m 


(1  -  u)  (flt)m+2  .  u  (at)11*2  ,  am+2  (t  +  l)"*2  ln( t  +  1) 
(m  +  l)(m  +  2)  (m  +  2) <*•  ^  (m  +  l)(m  +  2) 


_ Id _ 

(m  +  1)  (m  +  2) 


rrH-1 

Z 


k=l 


ak  (t  +  pk  (at)in~k+2 

(m  -  k  +  2) 


(15) 


Hysteresis 

The  hysteresis  strain  history  involves  a  constant  strain  rate  loading 
followed  by  an  unloading  at  the  same  rate. 


£(t)  =  Bt 

■  a ( 2tx  -  t) 


octet. 
tL  <  t  <  2Z. 


(16) 
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Substituting  into  the  constitutive  relation  yields  for  the  unloading 
history 


I 

I 


<« 


1 


V  • 


<r<t)  =  e  X  c 


m 


m 


r0(2t,-t)  m 

/  * 


-  f  [1  +  |i(l  +  t  -t)]  dt]  dr^ 


Cl  [i  +  \xa  +  t  -t)  dt 


(17) 


Integrating  yields 


0"(t)  =  E  C 


m 


m 


(1  -  -  fit) 


m+2 


pi[fl(2t,  -"t)]m+1 
(m  +  1) (m  +  2)  (m+  1) 


,m+2. 


i(2tl  "  fc)  in(l  +  2t  -  2t,)  M[8<1  +  fc)  3  tint1  +  t)  -  ln(l  +  2t  -  21^)] 

(m  +  2)  ^  +  ' 


(m 


,  ,  .»  T  [>a  +  tik  (26ti  -  at> 

4-  1)  (m  +  2)  j^L  („,  Util 


(m  +  1)  (m  +  2) 
m-k+2l 


+  (m  +  2) z  (2Btl  " 


fit) 


m-2 


-  1}  [6(2^  -  t)]™*1  0(1  +  t  -  t2)  ln(l  +  t  -  tx) 


(18) 


Relaxation 

Hie  relaxation  strain  history  may  be  written 


£(t)  =  0t 
=  atl 


0  <  t  <  t. 


The  constant  strain  rate  equation  is  valid  for  times  less  than  t^,  during 
the  strain  ranp.  For  the  relaxation  period,  t  >  t^. 


lo 


(T(t)  =  E 


m 


(flt,)”*2  (1  ~  _ 

t _ *  <  i  \  / _  i 


,mt2 


(m  +  1)  (m  +  2)  (m+l)(m+2) 


■fjr  [6(1  +  t)]“~  ln(l  +  t  -  t^) 


.  m+2 


m+2 


+  (m  +  1)  (m  +~2)  l(SV  md  +  t-V  +  (B{1  +  fc»  111(1  +  t)] 


u  V2  (Bt.)1*"**2^!  -I-  t)3k 

(m  +  l)  (m  +  2)  *  (m  -  k  +  2) 


-  (m  +  2) ^  (flV 


m+2 


t  > 


h 


(20) 


«4 

>■ 

a 


The  mathematical  model  was  used  to  predict  the  response  to  various 
strain  histories  and  compared  with  experimental  data.  The  distribution 
functions  for  the  human  spinal  anterior  and  posterior  longitudinal  liga¬ 
ments  and  the  supraspinous  ligament  are  shown  in  Figure  7.  These  functions 
were  determined  from  high  strain  rate  tests  by  numerical  solution  of  equa¬ 
tion  10.  These  distribution  functions  were  used  to  predict  responses  of 
the  different  ligaments  to  particular  strain  histories.  Carparisons  of 
theoretical  and  experimental  data  for  one  human  ligament  for  different 
strain  rate  histories  are  shewn  in  Figures  7  to  10. 
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Figure  3:  Failure  Stress  Versus  Number  of  Cycles  for  Thoracic  and  Lumbar  Ligamenta 
Flava  (L.F.)  from  Three  Rhesus  Monkeys  and  Two  Baboons.  The  L.F.  at  the 
thoracic  levels  (T3-T4,  T6-T7,  T9-T10)  and  lumbar  levels  (T12-L1,  L2-L3, 
L5-L6)  were  tested  in  rhesus  monkeys  A182 (Q),  N75 (A),  and  A322  (O)  and  in 
baboons  H40(H)  and  H54(#).  The  bar( — )  indicates  the  maximum  cyclic 
stress. 
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ALL-THORACIC  LEVEL 


Figure  5:  Failure  Stress  Versus  Number  of  Cycles  for  Thoracic  and  Lumbar  Anterior 

Longitudinal  Ligaments  (A.L.L.)  from  Three  Rhesus  Monkeys  and  Two  Baboons. 
The  A.L.L.  at  the  thoracic  levels  (T3-T4,  T6-T7,  T9-T10)  and  lumbar 
levels  (T12-L1,  L2-L3,  L5-L6)  were  tested  in  rhesus  monkeys  A182(Q), 
N75(A),  and  A322  (O)  and  in  baboons  H40®)  and  H54(#).  The  barf — )  indi¬ 
cates  the  maximum  cyclic  stress. 
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Figure  8.  100%/sec.  Strain  Rate  Response  of  a 

Human  Anterior  Longitudinal  Ligament  (Ll  -  L2) 
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Figure  9.  1%/sec.  Strain  Rate  Response  of  a 

Human  Anterior  Longitudinal  Ligament  (Ll  -  L2) 
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Table  2.  Single  Extension  of  Spinal  Ligaments  from  Rhesus  U70 


Cross- 


Spinal  Liqament 

Failure 

Load 

(N) 

Sectional 

Area 

(MIT) 

Failure 

Stress 

(MPa) 

Extension  at 
Failure 
(rnn) 

SSL 

C7/T1 

T3/T4 

9.1 

0.301 

0.394 

23.2 

2.24 

T6/T7 

15.6 

0.419 

37.1 

0.99 

T9/T10 

25.8 

0.477 

54.0 

1.04 

T12/L1 

11.3 

0.590 

19.2 

2.08 

L2/L3 

26.2 

0.247 

106.2 

2.44 

L5/L6 

15.6 

0.271 

57.5 

4.19 

LF 

C7/T1 

129.0 

6.72 

19.2 

3.76 

T3/T4 

83.2 

6.49 

12.8 

2.24 

T6/T7 

81.8 

4.94 

16.6 

3.66 

T9/T10 

242.0 

12.90 

18.8 

4.19 

T12/L1 

294.0 

14.90 

19.7 

3.68 

L2/L3 

702.0 

31.20 

22.5 

5.08 

L5/L6 

382.0 

25.30 

15.1 

5.36 

PLL 

C7/T1 

80.1 

2.41 

33.3 

0.76 

T3/T4 

71.0 

1.61 

44.0 

1.83 

T6/T7 

57.2 

0.90 

63.4 

1.12 

T9/T10 

T12/L1 

L2/L3 

19.6 

1.47 

1.96 

1.62 

11.5 

1.14 

L5/L6 

101. o 

1.31 

77.1 

1.68 

ALL 

C7/T1 

53.6 

4.83 

11.1 

1.75 

T3/T4 

136.1 

7.02 

19.4 

3.33 

T6/T7 

203.5 

10.3 

19.9 

3.12 

T9/T10 

494.6 

11.1 

44.7 

4.93 

T12/L1 

655.6 

10.8 

60.5 

3.53 

L2/L3 

381.6 

12.3 

30.9 

4.17 

L5/L6 

321.2 

10.1 

31.7 

3.79 

Table  3.  Single  Extension  of  Spinal  Ligaments  from  Rhesus  088 


Cross- 


Failure 

Sectional 

Failure 

Extension  at 

Load 

Area 

Stress 

Failure 

Spinal  Liqament 

(N) 

(imr) 

(MPa) 

(nm) 

SSL 

C7/T1 

32.9 

0.265 

124.5 

6.07 

T3/T4 

20.1 

0.367 

54.8 

5.46 

T6/T7 

36.5 

0.710 

51.4 

3.81 

T9/T10 

11.6 

1.280 

9.0 

3.30 

T12/L1 

16.0 

1.670 

9.6 

5.21 

L2/L3 

21.4 

0.678 

31.5 

4.45 

L5/L6 

37.5 

0.678 

55.3 

3.61 

LF 

C7/T1* 

2.64 

T3/T4* 

113.0 

6.25 

18.1 

5.84 

T6/T7* 

258.9 

11.40 

22.8 

4.06 

T9/T10* 

16.70 

T12/L1* 

489.3 

20.60 

23.8 

3.56 

L2/L3* 

658.3 

41.30 

15.9 

4.32 

L5/L6 

601.4 

41.30 

14.6 

4.29 

♦Sample  was 

split,  loads 

added,  extensions  averaged. 

PLL 

C7/T1 

43.3 

0.975 

44.5 

1.09 

T3/T4 

59.2 

1.480 

39.9 

1.78 

T6/T7 

1.460 

T9/T10 

1.540 

T12/L1 

24.9 

2.110 

11.8 

1.65 

L2/L3 

207.3 

0.301 

688.0 

1.52 

L5/L6 

114.8 

0.301 

380.9 

2.22 

ALL 

C7/T1 

90.5 

5.4 

16.9 

1.73 

T3/T4 

259.2 

8.6 

30.1 

2.41 

T6/T7 

386.1 

13.3 

29.0 

2.79 

T9/T10 

658.3 

18.5 

35.5 

3.94 

T12/L1 

581.8 

22.3 

26.1 

4.32 

L2/L3 

676.1 

13.3 

50.9 

4.19 

L5/L6 

770.4 

13.3 

58.0 

6.60 

Table  4.  Cyclic  Loading  of  Spinal  Ligaments  from  Rhesus  N75 


Cyclic  Cyclic 


Spinal 

Liganent 

Load  Ranqe 
Max  Min 
<N)  (N) 

Cross-Sectional 

Area 

(imr) 

Stress  Range 
Max  Min 

(MPa)  (**>a) 

Failure 

Stress 

(**>a) 

Number  of 
Cycles 

SSL 

T3/T4 

T6/T7 

0.37 

T9/T10 

0.32 

T12/L1 

35.4 

6.2 

0.27 

129.8 

22.7 

129.8 

3 

L2/L3 

0.45 

L5/L6 

49.4 

35.0 

0.52 

94.4 

66.9 

94.4 

2 

LF 

T3/T4 

T6/T7 

4.79 

T9/T10 

92.7 

46.3 

7.58 

12.2 

6.1 

12.2 

7,700 

T12/L1 

10.40 

L2/L3 

12.20 

L5/L6 

13.2 

PLL 

T3/T4 

T6/T7 

32.7 

23.7 

0.252 

129.8 

94.0 

129.8 

65 

T9/T10 

3.10 

T12/L1 

2.34 

L2/L3 

133.4 

41.9 

1.92 

69.6 

21.9 

95 

L5/L6 

1.65 

ALL 

T3/T4 

T6/T7 

208.2 

158.6 

8.18 

25.5 

19.4 

25.5 

3 

T9/T10 

281.8 

248.2 

11.40 

24.7 

21.8 

24.7 

343 

T12/L1 

423.9 

302.5 

11.80 

35.8 

25.6 

35.8 

3 

L2/L3 

5.74 

L5/L6 

390.1 

274.6 

4.90 

79.6 

56.0 

4,000 
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Table  5.  Cyclic  Loading  of  Spinal  Ligaments  frcm  Rhesus  A182 


Spinal 

Ligament 


Cyclic  Cyclic 

load  Range  Cross-Sectional  Stress  Range  Failure 
Max  Min  Area  Max  Min  Stress 

iHl  i*ii  (nnr)  (MPa)  (M>a)  (M>a) 


Number  of 
Cycles 


T3/T4 

0.437 

T6/T7 

15.3 

8.7 

0.559 

27.3 

15.5 

4,771 

T9/T10 

12.8 

3.3 

0.561 

22.8 

5.8 

9,100 

T12/L1 

24.1 

3.0 

0.561 

42.9 

5.4 

4,400 

L2/L3 

15.3 

2.5 

0.703 

21.7 

3.6 

22.1 

120 

L5/L6 

14.1 

2.5 

0.703 

20.1 

3.6 

29,160 

T3/T4 

3.8 

T6/T7 

4.6 

T9/T10 

11.2 

T12/L1 

389.9 

266.4 

16.7 

23.3 

15.9 

27.7 

1,100 

L2/L3 

378.8 

262.6 

23.6 

16.1 

11.1 

32,451 

L5/L6 

380.8 

260.1 

23.6 

16.1 

11.0 

19.8 

35,000 

T3/T4 

47.8 

32.9 

2.0 

23.6 

16.2 

5,700 

T6/T7 

43.3 

28.6 

2.0 

21.4 

12.7 

53.6 

4,300 

T9/T10 

2.8 

19.6 

o 

T12/L1 

3.3 

12.8 

0 

L2/L3 

2.2 

29.8 

0 

L5/L6 

66.5 

46.3 

2.2 

30.4 

21.1 

31.0 

300 

T3/T4 

6.8 

21.7 

o 

T6/T7 

198.7 

138.2 

13.8 

14.4 

10.0 

4,650 

T9/T10 

204.8 

146.6 

20.7 

9.9 

7.1 

17.5 

4,500 

T12/L1 

407.8 

289.6 

21.2 

19.2 

13.6 

10.000 

L2/L3 

436.0 

293.4 

13.1 

33.3 

22.4 

5,400 

L5/L6 

406.2 

174.7 

13.1 

31.0 

13.3 

36.0 

30,000 

Table  7.  Cyclic  Loading  of  Spinal  Ligaments  frcm  Baboon  H40 


Spinal 

Ligament 


Cyclic 
Load  Range 
Max  Min 

(N)  <N) 


Cross-Sectional 
Area 
(mm  ) 


Cyclic 
Stress  Range 
Max  Min 
(MPa)  (J*>a) 


Failure 

Stress 

(J*>a) 


Number  of 
les 


’4 

’7 

'10 

'Ll 

j3 

j6 

206.2 

287.3 

191.3 

1 

0 

1 

144.1  1 

211.9  2 

136.8  3 

’4 

’7 

143.9 

113.2  9 

’10 

JLZ 

9 

'Ll 

O 

482.6 

335.4  5 

1  0 

Jj 

j6 

661.9 

16 

486.5  26 

206.4 

112.9 

307.1 

220.3 

391.0 

287.3 

438.5 

303.5 

274.0 

202.1 

396.3 

287.7 

Table  8.  Cyclic  Loading  of  Spinal  Ligaments  frcm  Baboon  H54 


Cyclic  Cyclic 


Spinal 

Liqament 

Load  Ranqe 
Max  Min 
(N)  (N) 

Cross-Sectional 

Area 

(nnr) 

Stress  Ranqe 
Max  Min 

(MPa)  (MPa) 

Failure 

Stress 

(MPa) 

Number  of 
Cycles 

SSL 

T3/T4 

1.53 

T6/T7 

70.7 

52.6 

1.53 

46.2 

34.4 

46.2 

1520 

T9/T10 

69.8 

54.9 

2.53 

27.6 

21.7 

72.2 

4620 

T12/L1 

48.5 

16.3 

2.37 

20.5 

6.9 

34.8 

50 

L2/L3 

135.9 

103.8 

2.36 

57.6 

44.0 

63.9 

4057 

L5/L6 

68.9 

35.4 

2.19 

31.4 

16.1 

57.7 

26 

LF 

T3/T4 

10.90 

T6/T7 

10.10 

T9/T10 

17.20 

T12/L1 

301.3 

229.2 

25.20 

12.0 

9.1 

3994 

L2/L3 

51.80 

L5/L6 

72.70 

PLL 

T3/T4 

66.1 

48.5 

1.46 

45.4 

33.2 

49.7 

3982 

T6/T7 

66.6 

48.9 

1.73 

38.5 

28.3 

4000 

T9/T10 

36.3 

26.  L 

2.22 

16.4 

11.8 

92.4 

1240 

T12/L1 

2.43 

40.2 

0 

L2/L3 

70.7 

23.6 

9.04 

7.8 

2.6 

11.7 

22 

L5/L6 

195.3 

148.6 

8.42 

23.2 

17.6 

23.2 

258 

ALL 

T3/T4 

206.2 

159.5 

9.08 

22.7 

17.6 

22.7 

75 

T6/T7 

201.7 

159.5 

17.00 

11.9 

9.4 

28.1 

3985 

T9/T10 

207.5 

162.7 

24.80 

8.4 

6.6 

33.9 

3990 

T12/L1 

508.3 

396.8 

24.60 

20.6 

16.1 

36.9 

5743 

L2/L3 

525.9 

409.1 

30.20 

17.4 

13.6 

27.6 

3985 

L5/L6 

529.6 

414.5 

36.90 

14.4 

11.2 

34.0 

3985 

Table  9:  Smeary  of  Failure  Stress  for  Single  Extensions 
and  Cyclic  Loading  of  All  Animals. 


Combined 

Single  Extension  Cyclic  Loading  Conditions 

i  T  “ 1 t i  "  l 


Region 

Ligament 

Failure 

Stress 

(MPa) 

Var* 

(MPa) 

n 

Failure 

Stress 

(MPa) 

Var* 

(MPa) 

n 

Failure 

Stress 

(MPa) 

Var* 

(MPa) 

n 

Thoracic : 

S.S.L. 

50.6 

36.9 

7 

66.0 

79.7 

3 

55.2 

32.1 

10 

L.F. 

18.1 

3.3 

6 

13.4 

46.2 

14.5 

2 

16.9 

3.6 

8 

P.L.L. 

45.0 

11.2 

5 

49.2 

36.3 

10 

47.8 

29.2 

15 

A.L.L. 

25.8 

11.0 

8 

27.9 

17.2 

12 

27.1 

14.8 

20 

Lumbar: 

S.S.L. 

46.6 

34.9 

6 

71.9 

43.0 

12 

63.5 

41.2 

18 

L.F. 

18.6 

4.0 

6 

31.6 

24.7 

7 

25.6 

18.9 

13 

P.L.L.** 

[233.9] [293.6] 

5 

54.0 

34.0 

12 

[106.9] [172.9] 

17 

A.  L.  L. 

43.0 

15.2 

6 

37.9 

12.8 

11 

39.7 

13.4 

17 

♦Variability  (See  text). 

♦♦Brackets  indicate  questionable  stress  values  based  on  small  and  difficult  to 
delineate  ligament  areas. 
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Table  10:  Cyclic  Extension  at  Failure 


Combined 

Single  Extension  Cyclic  Loading  Conditions 

Failure  Failure  Failure 


Region 

Ligament 

Extensn 

(mm) 

Var* 

(itm) 

n 

Extensn 

(mm) 

Var* 

(mm) 

n 

Extensn 

(mm) 

Var* 

(mm) 

n 

Thoracic : 

S.S.L. 

3.27 

2.01 

7 

2.43 

3.28 

3 

3.02 

1.73 

10 

L.F. 

3.96 

1.16 

6 

3.26 

1.81 

1 

3.86 

1.09 

7 

P.L.L. 

1.32 

0.47 

5 

1.34 

0.80 

9 

1.33 

0.68 

14 

A.L.  L ■ 

3.00 

1.09 

8 

2.50 

0.62 

7 

2.76 

0.91 

15 

Lumbar: 

S.S.L. 

3.66 

1.21 

6 

3.14 

1.62 

7 

3.38 

1.41 

13 

L.F. 

3.88 

1.69 

6 

4.88 

1.88 

5 

4.34 

1.77 

11 

P.L.L. 

1.64 

0.39 

5 

2.00 

1.16 

10 

1.88 

0.97 

15 

A.L.  L. 

3.77 

2.02 

6 

2.94 

1.34 

9 

3.27 

1.63 

15 

♦Variability  (See  text) 
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